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I. INTRODUCTION 

There exist two well established methods to trap charged particles: the Penning trap [1] and the Paul trap [2]. In the Penning 
trap the particle is confined in space by a combination of static magnetic and electric fields. In the Paul trap the trapping is 
caused by a high frequency electric quadrupole field. The subject of this article is to present a third mechanism for trapping 
charged particles — trapping by beams of electromagnetic radiation. It was discovered some time ago [3-6] that a properly 
prepared beam acts as a "waveguide" for particles, confining their motion in the transverse directions. Similar phenomena occur 
in RF fields [7]. In all these cases, the essential role is played by the electric field configuration in the plane perpendicular to 
the beam axis (for nonrelativistic electrons, the magnetic field is less important). Particles are confined to the vicinity of the 
minimum-energy points. In particular, for beams of electromagnetic radiation carrying orbital angular momentum such points 
lie on the beam axis. One beam may confine particles only in the transverse direction. Two or three crossing beams may fully 
confine particles, acting as a substitute for the Paul trap. 

Trapping of charged particles by beams of electromagnetic radiation is based on the same general principles as the Paul trap. 
In the Paul trap the quadrupole radio-frequency field is produced by a system of electrodes. In our case, the field is that of freely 
propagating electromagnetic beams. In both cases the essential role is played by two factors: the proper shape of the force field 
and fast oscillation or rotation. These features are clearly seen when we invoke the notion of the ponderomotive potential. This 
notion is applicable when the the motion of a charged particle can be unambiguously split into a slow drift and a fast quiver 
motion. Such a splitting is possible when a nonrelativistic particle is moving in a rapidly oscillating field. We shall restrict 
ourselves in this paper to monochromatic beams since they are especially well suited for trapping particles. The best example of 
such beams are Bessel beams which will be analyzed in detail. 

II. ELECTRIC AND MAGNETIC FIELDS OF BESSEL BEAMS 

Every monochromatic electromagnetic field can be represented by the following general expressions: 

E(r,t) = E c (r) cos(ujt) + E s (r) sin(cjt), (la) 
B(r,t) = B c (r)cos(u>t) + B s (r)sm(wt). (lb) 

According to Maxwell equations the electric and magnetic field vectors are interdependent, 

u:B c {r) = V x E s (r), uB s (r) = -V x E c (r). (2) 

Generically, the electric and magnetic field vectors at each point trace two ellipses spanned by the pairs of vectors {E c , E s } and 
{B c , B s }. 

Bessel beams are most easily analyzed with the use of the Riemann-Silberstein (RS) vector [8], defined as a complex combi- 
nation of the electric and magnetic fields, 

F(r,t)=E(r,t)+icB(r,t). (3) 

This vector obeys the following equations: 

d t F{r,t) = cV x F{r,t), V-F(r,i)=0. (4) 

For monochromatic fields we can write 

F(r,t) = e - ,ut F + (r) + e mt F-(r), (5) 

and the Maxwell equations reduce to the two equations for Trkalian fields [9] 

cV x F ± (r) = ±ujF±(r). (6) 

The two parts in (5), as we show below, correspond to two helicities. By comparing the formulas (la) and (5), we find 

E c = RcF + + ReF_ , E s = lmF + - ImF_ , (7) 
B c = RcF + + ImF_ , B s = -ImF + + RcF_ . (8) 

It follows from these formulas that for fields of definite helicity the ellipses traced by the electric and magnetic field vector 
coincide. 
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Following Ref. [10] (with some small changes in notation), we shall describe the electromagnetic field of the Bessel beam in 
cylindrical coordinates, with the z axis in the direction of propagation, in terms of the RS vector in the following form: 



-ixut & ik\\z 



F x (p,<p,z,t) 
F y (p,(f>,z,t) | = E e 
F z {p, <j),z,t) 

K _ e i(M + D4>j M+l{k±p) + K+e m-i^j M _ l{k±p) 

x | -i K _e i ( M + 1 ^J M +i{k±p)+iK + e^ M - 1 ^J M -i(k±p) 
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(9) 



where 



x Jkl + kl±k ll 



k± 



(10) 



We shall consider only positive values of M since Bessel beams with negative values can be obtained (up to an overall sign) by 
the following transformations of spacetime variables and helicity: 



M 



-M, x^-X. 



-(f), t ->■ -t. 



(11a) 



The formulation in terms of the RS vector enables one to introduce, by analogy with quantum mechanics, the classification of the 
solutions of Maxwell equations with the use of quantum numbers [10]. One may view these quantum numbers as characteristics 
of the photons making the beam. In the case of Bessel beams these quantum numbers are the eigenvalues of the following 
four operators (we dropped h in the definitions): the component of momentum in the direction of propagation k», the length of 
momentum in the transverse direction k±, the projection of the total angular momentum (orbital plus spin) on the direction of 
propagation M z , and the helicity operator (the projection of the total angular momentum on the direction of momentum) x, 



= -id z , k± = \J-dl - d$, 



(12) 



M z = -i(xd y - yd x ) + 



-i 
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(13) 
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-d x > dy, 

d z > -d x . 
-d y ' d x > 



(14) 



The 3x3 matrix in Eq. ( 1 3) represents the z component of the photon spin in the vector representation. The representation of the 
helicity operator in the form (14) was given in [1 1]. The factors exp[z(M± l)(f>] and cxp[iM<f>] in the definition of Bessel beams 
reflect the decomposition of the total angular momentum M into the components with three different orbital angular momenta. 

In classical physics quantum numbers are just convenient labels attached to the solutions of Maxwell equations. Thus, each 
Bessel beam (9) is characterized by the following four parameters: a real parameter fcii, a positive parameter k±, an integer M, 

and x (±1)- Every Bessel beam is monochromatic with the frequency uj = c^Jk^ + k 2 ± . The amplitude E in Eq. (9) measures 

the strength of the beam in volt per meter (1 V/m corresponds to the beam intensity 2.65 x 10 -3 W/m 2 ). The dimensionless 
electric and magnetic field vectors for the Bessel beam are obtained by taking the real and imaginary parts in Eq. (9) and dropping 
the factor E n , 
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FIG. 1. Electric field vectors in the xy plane near the beam axis for \ = 1 an d f° r Af = 0, 1, 2, 3. This field configuration corresponds to 
z — and t = 0. For other values the vectors would be rotated by the angle (fcy z — x^i). The top configurations correspond to M = and 
M — 1. The bottom configurations correspond to M = 2 and M = 3. 



where 

c M = cos(fc||z - x^t + M4>)J M (k±p) , (16a) 
s M = sin(fc||Z - x^t + M<fi)J M (kj_p) . (16b) 

In Fig. 1 we show the electric field configuration in the transverse plane near the beam axis for different values of M. Note that 
for M = 2 the configuration of electric forces is the same as of the gravitational forces on the saddle surface in the mechanical 
model of the Paul trap [2]. Therefore, at least in this case we can expect trapping of particles in the transverse plane. To exhibit 
rotation of field vectors in time, we plotted these vectors in Fig. 2 at the time one quarter of the cycle later. In Fig. 3 we show 
the rotation of the field vectors for M = 2. 



III. MOTION OF CHARGED PARTICLES IN A BESSEL BEAM 



All field components of the Bessel beam with \M\ > 1 vanish along the z axis since in these cases all Bessel functions in 
Eq. (9) vanish at the origin. The motion of charged particles in Bessel beams in the vicinity of this line was thoroughly studied 
in Ref. [4] by solving numerically the standard equations of motion 

m^=e/-(*( S ))%^, (17) 
as z as 

where s is the proper time. These equations can be rewritten in the following dimensionless form 

-W_ 7 ^ ({(r)) *M 

where the particle coordinates £ M are measured in units of the longitudinal wave length 2n/ku and the dimensionless proper time 
t is measured in units of the inverse circular frequency t = ujs. The dimensionless effective coupling constant 7 is 

7 _ eg ° fc H _ eE *l k W (19) 
raw 2 mc 2 (l + (fcj_/fc||) 2 ) 

The dimensionless electromagnetic field tensor ip^ is composed of the dimensionless electric and magnetic components given 
by Eqs. (15). The parameter 7 is usually very small. It is approximately equal to the ratio of the energy gain (i.e. the energy that 



5 






1 


\ 






1 


\ 








\ 


. ~ ^ s 


\ X - - 
\ \ « ■ 






- - , / 
■'ft 


t t ' ■ 






. « \ \ 








- - X \ 


s ^ ~ - 










\ 


1 






\ 


{ 





FIG. 2. Electric field vectors in the xy plane as in Fig. 1 but a quarter of the cycle later. Each field vector is rotated counterclockwise by ty/2. 
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FIG. 3. Snapshots of the electric field vector in the transverse plane for M = 2 and x = 1, taken at t = 0, T/4, T/2, 3T/4, where T is the 
period of the wave. 



the particle acquires moving under the influence of the electric field Eq along the distance of 1/fcy) to the particle rest energy 
mc 2 . The particle trajectories scaled by the wavelength with initial velocities scaled by the speed of light will look the same 
for different situations provided the value of 7 remains the same. The smaller the mass, the more effective is the action of the 
electromagnetic wave. Therefore, in all our specific examples we shall consider electrons. In what follows we shall consider 
nonrelativistic electrons whose motion is governed mainly by the electric field. One beam may confine the motion of an electron 
only in the transverse directions — the electron is guided along the direction of propagation. For full trapping we will need at 
least two beams. 

Typical trajectories of an electron guided by the Bessel beam are shown in Fig. 4. For Bessel beams describing the propagation 



FIG. 4. Three trajectories of an electron guided by the Bessel beam with M — 2 and \ — 1 obtained for different initial transverse velocities. 
All distances are measured in units of the longitudinal wavelength. Thick line represents the beam axis. Note a large difference in the scaling 
of transverse and longitudinal distances. 




FIG. 5. Four trajectories of an electron guided by the Bessel beams with M = 0, 1, 2 and 3. The tightly trapped trajectories for M = and 
M = 2 are indistinguishable. Trapping is less pronounced for M = 3 (wide trajectory). There is no trapping for M = 1 (runaway trajectory). 



of light, the effective transverse dimension 2ir/k± is much larger than the longitudinal wave length. This means that the 
dimensionless parameter fcj_/k|| is usually very small. The trajectories shown in Fig. 4 are obtained for k±/k\\ = 0.01 and 

7 = 0.00001. For visible light, this value of 7 corresponds to the beam intensity of the order of 10 13 W/m 2 . This high intensity 
resulted from our choice of rather high transverse velocities that were taken for illustrative purposes. The three trajectories 
guided by the beam start at the same point near the beam axis and with the same initial longitudinal velocity z(0) = 0.001c 
but they have different initial transverse velocities x(0) = 0.0002c, 0.0001c, and 0.00005c. Trapping of electrons with initial 
transverse velocity exceeding the value 0.001c will require the increase of the beam intensity. For moderate beam intensities, 
relativistic electrons will not be trapped. 



FIG. 6. Projections on the transverse plane of the three trajectories of an electron guided by the Bessel beams with M = 2. The trajectories 
started with the same initial velocities in the vicinity of the zeros of the Bessel function. 



In Fig. 5 we plotted the trajectories trapped near the beam center for different M obtained when x(0) = 0.0001c and with 
the values of the remaining parameters the same as in Fig. 4. We can easily explain qualitatively these results, with the help of 
the formulas (15) and Fig. 1. First of all, when k±/ku is small, only the Bessel function Jm- x contributes significantly to the 
transverse components of the fields. This observation underscores the role of the orbital angular momentum. For realistic Bessel 
beams only one orbital angular momentum component, with M — x, plays the dominant role. Therefore, for \ = 1, the trapping 
will occur only for M — 0, 2 and 3 but not for M = 1. This is so because in the last case the dominant part comes from the 
Bessel function Jo an d near the beam center the field configuration is that of a circularly polarized plane wave. Trapping is the 
best for M — and M — 2 and for \ = 1 because only in these two cases the dominant role in the formula (9) is played by 
the Bessel function J\. Therefore, the field strength increases linearly with the distance from the beam axis, as in the Paul's 
mechanical model. For M — 3 the trapping is weak, because the Bessel function J 2 increases more slowly near the beam center, 
so that the radius of the effective potential well is larger. With the increase of M, the trapping will become weaker and weaker. 

The problem of an electron motion near the Bessel beam axis was solved explicitly in Ref. [3] in classical and in quantum 
mechanics for the fields described in the paraxial approximation. The analytic solution in this case was possible due to the fact 
that the field configuration close to the center of the beam with M = 2 and x = 1 is particularly simple — the field increases 
linearly with the distance from the beam axis. The main feature of the analytic solutions is the separation of the longitudinal 
motion from the motion in the transverse plane. The motion in the transverse plane is bounded whereas the longitudinal motion 
decomposes into oscillations and a uniform motion with velocity determined by the initial conditions. The analytic solutions of 
the equations of motion (17) in this simplified model fully confirm our numerical results. 

Bessel functions vanish not only at the center of the beam but they also have a succession of other zeros. It turns out that 
charged particles may be trapped in regions corresponding to these zeroes. In Fig. 6 we show the four trajectories trapped at the 
first four zeroes for M — 2 and \ = 1. 



The trajectories of particles trapped by Bessel beams in the nonrelativistic regime can be reproduced remarkably well in 
the ponderomotive approximation. The notion of the ponderomotive potential introduced by Gaponov and Miller [12] and 
independently by Boot and Harvie [13] is very useful when the particle moves slowly in an inhomogeneous, rapidly oscillating 
electromagnetic field. The approximation is based on the decomposition of the particle trajectory into a sum R{t) + r(t), where 
the first part is slowly varying and the second part describes small but fast oscillations with zero mean value. We start from the 
equations of motion with the Lorentz force, 



The necessary condition for the applicability of the ponderomotive approximation is \v ■ V/|/|w/| <C 1, where / is the force 
field. Under this assumption, one may expand the fields in the equations of motion around the slow trajectory R(t) keeping only 
linear terms in r{t). In the lowest order, using the notation of Eqs. (la) and (2), we obtain 



IV. PONDEROMOTIVE POTENTIAL 




(20) 



r = — [E C (R) cos(u;i) + E S (R) sin(wi)] . 



(21) 
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After disregarding the slow variation of R with time, we obtain by a straightforward integration 

e 



[E C (R) cos(wt) + E S (R) sin(wt)] 



(22) 



where we dropped the integration constants to obtain zero mean values. In the next order, the equation reads 



R 



— \\(r ■ V)E C (R) + -r x (V x E S (R)) 1 cos(o;t) 
m U uj . 

[(r • V)E S (R) - -j-r x (V x £7 C (JR)) sin(urt)}. 



(23) 



In the final step, after inserting the solutions of Eq. (21) for r and r, we perform time averaging of the right hand side over one 
period of fast oscillations 2it/uj and we obtain 

R - -7r4-T \ ( E c{R) ■ V) E C (R) + E S (R) x (V x E S (R)) 



+ (E S (R) ■ V) E S (R) + E C (R) x (V x E C (R)) 

9 

V [E c (Rf + E S (R) 2 ] = -VV P (R), 



Am 2 oj 2 



where we used of the vector identity 



E x (V x E) = -VE 2 -(E-V) E. 



The ponderomotive potential V P (R) can also be expressed as the square of the electric field averaged over time 



V P (R) = 



(E(R,t)-E(R,t)). 



(24) 



(25) 



(26) 



Despite the fact that only the electric field enters into this formula, the ponderomotive potential contains also a contribution from 
the magnetic field, as seen from our derivation. In our dimensionless variables the equations of motion have the form 



7 



V<5«,r).£(€,r)>. 



(27) 



The use of ponderomotive approximation requires some care in handling the initial conditions. There is no need to adjust the 
initial value of the slowly varying part of the trajectory (the term R in Eq. (21)). We can always identify R(t = 0) with the 
initial position of the exact trajectory because the amplitude of fast oscillations is small. However, the adjustment of the initial 
velocity is necessary because the velocity of fast oscillations is large and it gives a significant contribution to the true velocity. 
This adjustment is especially important when the trajectory starts in a region where the field is strong. According to Eq. (22), r 
can be related to the electric field by the following formula 

[E C (R) sm(ujt) - E S (R) cos(wi)] 



TOW 



; [E C (R) cos(wi + tt/2) + E S (R) sm(u>t + tt/2)] 
;E(R,t + Tr/2/u). 



(28) 



Since the ponderomotive approximation reproduces only the slowly changing part of the trajectory, the initial velocity of the slow 
motion R must be identified with the difference of the exact solution and r. The significance of the initial velocity adjustment 
will be illustrated with the examples given in the next section. 

The ponderomotive potential for a Bessel beam is a function of the distance from the axis only, 



v p(p) = v o (n + J M -i(kj_p)) 2 + 2 {J M {k±p)f + («;_ J M +i(fc_L/o)) 2 



(29) 



where V is the amplitude determined by the beam intensity; in our dimensionless variables, V = 7 2 /2. For realistic Bessel 
beams, when k± -C fcii, only one term in this formula contributes significantly. This term represents the contribution from the 
component with orbital angular momentum M — x- The ponderomotive potentials for \ = 1; A|| = 0.01, and for different 
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FIG. 8. Ponderomotive potential for the Bessel beam with M = 2 and \ — 1 as a function of x and y. 



values of M are shown in Fig. 7. The minima correspond to the zeros of the dominant Bessel function. Only the minima at 
p = for M > 1 are exactly equal zero, the remaining ones are close to zero. 

The trajectories obtained by integration of exact equations, depicted in Fig. 4 and in Fig. 5 coincide with those obtained from 
the ponderomotive potential (29). The vertical line in those figures is the line of the potential minima (except for M = 1 where 
it is a maximum). For M > 1 the beam axis is an obvious potential minimum, since the whole electromagnetic field vanishes 
there. For M — the z-component of the field does not vanish on the axis but its contribution is very small when k±/kn is small. 
The notion of the ponderomotive potential enables one to understand the existence of a limiting transverse particle velocity that 
allows for trapping. The transverse kinetic energy cannot exceed the height of the ponderomotive potential. 

The ponderomotive potential for the Bessel beams, in addition to the potential well along the axis, has a succession of potential 
valleys corresponding to the minima in Fig. 7. In Fig. 8 we show the potential for M — 2 and \ — 1 as a function of x and 
y where the potential valleys are clearly seen. The trajectories depicted in Fig. 6 are trapped in those valleys and they are 
reproduced very well by solving the equations with the ponderomotive potential. 

So far we considered trapping by beams carrying definite angular momentum. Such beams are most natural in this context 
since the ponderomotive potential has easily identified minima. However, as shown in the next section, trapping is found also 
for field configurations that do not have definite angular momentum. 



10 




FIG. 9. Contour lines of the ponderomotive potential in the xy plane created by a superposition of the Bessel beam with M — 2 and a plane 
wave with the same frequency lu drawn for four values of z. The dashed circle represents the projection of the helix on the xy plane. 




FIG. 10. Motion of a slow particle guided by the helix (thick line) created by a superposition of the Bessel beam and a plane wave. There are 
no noticeable differences between the exact solution (left) and the ponderomotive approximation (right). 



V. TRAPPING OF PARTICLES BY SUPERPOSITIONS OF BESSEL BEAMS 



A. Trapping of particles along a helix 

We shall now consider the motion of an electron in the superposition of a Bessel beam and a circularly polarized plane wave. 
We assume that both waves propagate in the same direction and have the same frequency. The dimensionless electric and 
magnetic fields are in this case 

acos(w(z/c — t) + (k-Cm+i + k+Cm-i) 

-asin(cj(z/c - t) + (k-s m +i - k+Sm-i) ), (30a) 
2sm 

asin(w(z/c - t) + (k_s m +i + k+Sm-i) 

acos(ui(z/c - t) + (-K-Cm+i + k+cm-i) ), (30b) 
-2c M 

where a measures the relative strength of the plane wave as compared to the Bessel beam. The ponderomotive potential is given 
by the formula 

V p (p, <j), z) = Vo[ (k+ JM-i(k±p)) 2 + 2 (JM(k±p)) + (n-JM+i(k±p)) 

+ a 2 + 2aK+J M -i(k±p) cos ((M - 1)0 - (w - k\\)z) } . (31) 
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FIG. 11. Motion of a fast particle in the field as in Fig. 10. There is a clear difference between the exact (left) and the approximate solution 
(right) but the ponderomotive potential still reproduces qualitatively the trapping near the helix. 




FIG. 12. Contour lines in the xy, yz, and zx planes of the ponderomotive potential for two crossed Bessel beams with M = propagating in 
the x and y directions. There is a potential minimum at the crossing point. 



The minima of V p (p, (f>, z) can be easily found. The minimum in variables <j> and z is reached when the cosine attains the lowest 
value -1. The minimum in p for this value of the cosine can be easily found numerically. These condition define the helix 
z = (u> — fcii ) ((f) + n)/ (M — 1). The line of the helix traces the minima of the ponderomotive potential and is the analog of the 
beam axis for a Bessel beam. The helical motion of these minima along the beam axis is visualized in Fig. 9. Therefore, in this 
case we can expect that the electrons will be guided along the helical line. However, this requires stronger fields. For electrons 
with initial transverse velocities the same as for a single Bessel beam we increased 7 to the value 0.0005. In Fig. 10 we show 
the trajectories for M — 2, \ = 1> an d a — 3 obtained by solving exact equations and those with the ponderomotive potential. 
The small parameter fcj_/fc|| determines the radius of the helix and to make this figure more transparent we have increased this 
parameter tenfold from the value used in previous figures. 

In order to trap more energetic particles, we need stronger fields. As an example, we have chosen the initial transverse velocity 
as 0.01c and we have increased the field strength to 7 = 0.01. In Fig. 1 1 we show that the particle is trapped around the helix. 
The ponderomotive approximation is still applicable but it does not reproduce exactly the fine details of the trajectories. 
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FIG. 13. Contour lines in the xy, yz, and zx planes of the ponderomotive potential for two crossed Bessel beams with M = 2 propagating in 
the x and y directions. There is a potential minimum at the crossing point. 




FIG. 14. Trajectories of a particle trapped in the vicinity of the potential minimum. Exact trajectory (left) is completely different from the 
trajectory obtained in the ponderomotive approximation without the velocity adjustment (right). 



B. Trapping of particles by crossed Bessel beams 



In all previous cases the trapping occurred only in the transverse direction. The trapping in all directions can be achieved 
by using crossed Bessel beams. Crossed Bessel beams trap particles in the vicinity of the crossing point. We have chosen two 




FIG. 15. The trajectory obtained in the ponderomotive approximation with the velocity adjustment (right) is much closer to the exact trajectory 
(left) than in Fig. 14. 
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FIG. 16. Contour lines in the xz plane of the ponderomotive potential for the standing Bessel wave (M = 0). Owing to cylindrical symmetry 
the contour lines would be the same after an arbitrary rotation around the beam axis shown as the center line. The minima of the potential are 
marked by the dots. 



identical Bessel beams propagating in the x and y direction. In this case, the ponderomotive potential is given by the formula: 

V p = V [k 2 + (Jm_i(Px) + Jm-i(Pv)) + K- (Jm+i(Px) + Jll+liPy)) 

- 2k + J m (Px)Jm-i(p v ) cos(ip + <j>y) + 2k-J m {Px)Jm+i{p v ) cos(ip - <f> y ) 

- 2K + J M -i(px)JM(p v )sm(xp - <j> x ) - 2k-J m +i(Px)Jm(p v ) sm(ip + <f> x ) 

- k\Jm-i(Px)Jm-i(p v ) sia(ip - <j> x + <j> y ) 
+ K 2 _J M+1 (p x )J M +i(p y )sin(ip + (j) x - (j>y) 

- JM-l(Px)JM+l(Py) Sm(lp - (f> x - 4>y) 
+ JM+l(Px)JM-l(p y )sm(xp + 4> x + 4>y) 

+ 2(j 2 M (p x ) + J 2 M (p y ))], (32) 
where we used the following abbreviations: 



<t> x = arctan(z/y), 
Px = k± \Jy 2 + z 2 , 



<p y = arctan(x/z), 
p y = kj_ \/ z 1 + x 2 . 



ip = fc|| (x-y) + M((j> x - <p y ) 



(33) 



Again the value M = 1 is exceptional — it does not lead to trapping. The shapes of the potential wells for M = and M = 2 
are very similar, as shown in Fig. 12 and Fig. 13). These shapes clearly indicate that the potentials produced by crossed Bessel 
beams can trap particles in the vicinity of the beam crossing point. Solutions of the equations of motion fully confirm these 
predictions. For higher values of M trapping is less effective since the wells become shallower. 

For crossed Bessel beams the ponderomotive approximation is less successful than for a single Bessel beam. In Fig. 14 and 
Fig. 15 the exact trajectory for M = is compared with the approximate trajectories. Even though the adjustment of the initial 
velocity significantly improves the result, there remain substantial differences. The discrepancy is most likely due to a rather 
complicated shape of the potential well which may lead to a chaotic behavior — small changes produce large effects. 



C. Trapping of particles by a standing Bessel wave 



In all previous cases we considered trapping by running waves. In practical applications it may be easier to work with standing 
waves. It turns out that a standing Bessel wave can also fully trap charged particles. A standing Bessel wave is constructed by 
adding two Bessel beams that differ only in the sign of fc||. The ponderomotive potential of this superposition is given by the 
formula 

V p (r) - 2V [2k 2 /k 2 ± (Jl_ 1 (k ± p) + J 2 M+1 (k±p)) 

+ (1 + cos(2fc||z)) (Jj f _ 1 {k±p) + J 2 M +i(k±p) + 2J 2 M (k ± p)) ] . (34) 
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FIG. 17. Trajectory of a particle trapped in the vicinity of the potential minimum at z = it/2k\\. Thick line is the beam axis. There is no 
noticeable difference between the exact solution (left) and the one obtained within the ponderomotive approximation (right). 




FIG. 18. Trajectories of a particle trapped in the vicinity of the potential minimum with the same parameters as in Fig. 17 but with the starting 
point moved away from the minimum. The trajectories are completely different since the adjustment of the initial velocity has not been made. 




FIG. 19. Trajectories of a particle trapped in the vicinity of the potential minimum with the same parameters as in Fig. 17 but with the starting 
point moved away from the minimum. The agreement between the exact trajectory and the one in the ponderomotive approximations has been 
fully restored by adjusting the initial velocity. 
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The ponderomotive potential is a periodic function of z with the period equal to 7r/fc||. The amplitude of these oscillations is 
exceedingly small. Owing to the large value of k^/k\, the ^-independent terms dominate. The case M — is exceptional since 
then the z-dependent term dominates near the beam axis. In Fig. 16 we show contour lines of the ponderomotive potential for 
M = cut by the xz plane that exhibit the potential minima. The trajectories of a particle trapped in three dimensions near 
such a minimum are shown in Figs. 17-19. In Fig. 17 we have a fairly good agreement between the exact trajectory and its 
ponderomotive approximation because the starting point is very close to the potential minimum. In Fig. 18 the trajectories are 
completely different because the adjustment of the initial velocity has not been made. Finally, in Fig. 19 the adjustment has been 
made and the agreement has been fully restored. 

VI. OUTLOOK 

We have shown that Bessel beams might be used as a "waveguide" and also as a trap for charged particles. So far we 
only considered a single charged particle. The extension of our results to many particles is an open problem because it would 
require the inclusion of the Coulomb interaction between particles. The standing Bessel wave configuration offers a possibility 
of trapping many particles in successive minima of the ponderomotive potential. In this way one may construct an analog of 
the optical lattice for charged particles. The only condition for such a construction is that the Coulomb repulsion should not 
overcome the ponderomotive potential barrier. This means that the intensity of the beam must satisfy the inequality I > I c = 
4:irm(c/ A) 3 . For the wave length A = 6 • 10~ 7 m, this gives a rather high (but not totally unrealistic) value, I c = 1.4 ■ 10 15 W/m 2 . 
When the intensity obeys the condition I < 2I C then only one electron can be trapped at each potential minimum — the Coulomb 
repulsion will push extra electrons away. One may significantly improve our standing wave trap by adding a constant magnetic 
field in the direction perpendicular to the beam axis. Such a magnetic field will also quantize the spin degree of freedom of 
electrons, transforming each electron into a qubit. All these conditions apply to electrons at rest. In order to hold electrons that 
undergo motion in a trap, we need higher intensities. 
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